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ABSTRACT 

We review the recently proposed light-cone quantization of the matrix model 
which is expected to have a critical point describing 2-d quantum gravity coupled 
to c = 2 matter. In the N — > oo limit, we derive a linear Schroedinger equation 
for the free string spectrum. Numerical study of this equation suggests that the 
spectrum is tachyonic, and that the string tension diverges at the critical point. 



★ Based on the talk delivered by I. R. Klebanov at the 7th Nishinomiya- Yukawa Memorial 
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1. Introduction 



Recently there has been considerable renewed interest in large-iV matrix models. The 
one-dimensional hermitian matrix chain models have been solved in the double-scaling limit 
and identified with the c < 1 minimal models coupled to two-dimensional quantum gravity [1, 
2] . In the same fashion the hermitian matrix quantum mechanics, which describes the c = 1 
theory, has been solved both for non-compact [3] and for circular target space [4] . The c = 1 
model can be interpreted in terms of D = 2 string theory [5], where the role of the extra 
dimension is played by the conformal factor of the world sheet quantum gravity. One should 
keep in mind, however, that the models with c > 1 are of much greater interest because they 
are expected to correspond to string theories in D > 2. These theories should have a much 
richer structure than in D = 2 because strings can exhibit transverse oscillations. In spite 
of some interesting new developments [6, 7], there is little solid information about the c > 1 
theories. 

From the existing analytical and Monte Carlo studies there is some evidence that the 
discretized c > 1 models are in the branched polymer phase [8] and do not lead to acceptable 
bosonic string theories. On the other hand, the continuum approach indicates that these 
models are tachyonic [9]. It is often suggested that the branch polymer behavior is indicative 
of the presence of the tachyon. However, some of the discretized models should not be 
tachyonic due to reflection positivity [8]. We believe that a deeper understanding of the c > 1 
models is needed because, even though the simplest such theories are probably unacceptable, 
there may be modifications that lead to new interesting string models. With this idea in 
mind, Dalley and one of the authors [10] have recently attempted a new approach to the 
theory of random surfaces embedded in two dimensions (c = 2). The corresponding matrix 
model is a two-dimensional super-renormalizable scalar field theory which can be quantized 
using the light-cone coordinates. The iV — > oo limit naturally leads to a free string light-cone 
Schroedinger equation [11]. In these notes we will review the light-cone quantization of this 
model, and will report on new numerical studies of the free string spectrum which suggest a 
different kind of behavior than what was anticipated in ref. 10. 
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2. Random Surfaces Embedded in One and Two Dimensions 



Before proceeding to the c = 2 model, we will briefly review some of the lessons learned 
from the exactly soluble c = 1 matrix model for the purpose of comparison. The discretized 
approach to 2-d quantum gravity coupled to c = 1 matter [12] is generated by the euclidean 
matrix quantum mechanics with the action 

J \2\dx J 2a' 3VN J V ; 

where M(x) is an N x N hermitian matrix. The connection of this matrix model with 
triangulated random surfaces follows, as usual, after identifying the Feynman graphs with 
the graphs dual to triangulations. The lattice link factor is the one-dimensional scalar 
propagator, 

G(xi,xj) = exp (—\x{ - Xj\/Va^ (2) 

so that the parameter a' sets the scale in the embedding dimension. The model (1) possesses 
a global SU(N) symmetry under M — > WMQ. It is well-known that the SU(N) singlet 
spectrum is described by N free fermions moving in the cubic potential, with the Planck 
constant ~ 1/N. In the leading WKB approximation the ground state energy is ~ iV 2 , 
which corresponds to the sum over surfaces of spherical topology. This sum is singular as 
A ~ A c — A — > 0, which is the value of the coupling where the Fermi level reaches the local 
maximum of the cubic potential, 

1 iV 2 A 2 

Z = 



/a 1 hiA 

The geometrical meaning of this singularity is that the Feynman graphs with very large 
numbers of vertices are becoming important. Near this singularity the universal continuum 
limit of the 2-d quantum gravity coupled to c = 1 matter can be defined. In this limit, the 
gaps in the low- lying spectrum go to zero as 1/| In A|. The underlying dispersion relation [5] 

a'E 2 -pl = 

is characteristic of a massless two-dimensional field theory. The hidden 0-dimension derives 
from the fluctuating conformal factor of the world sheet. It is worth emphasizing that 
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the string scale a! does not undergo any singular renormalization in the A — > limit and 
essentially corresponds to the "bare" string scale defined by the link factor (2). 

An important fact about the c = 1 matrix model is the decoupling of the SU (N) non- 
singlet states, whose energies diverge as | In A | . As a result, only the singlet states are 
important for random surfaces embedded in a non-compact dimension, or in a circle of 
radius R > R c [4]. The decoupling of the non-singlet states corresponds physically to the 
confinement of Kosterlitz-Thouless vortices [13, 14]. 

In ref. 10 a generalization of the matrix model (1) to two dimensions was studied, 

S E = J <PxTr Q(0«M) 2 + i/iM 2 - ^AM 3 ) , (3) 

where M(x°, x 1 ) is an N x N hermitian matrix field. Just as for c = 1, the Feynman graphs 
of this theory are identified with the graphs dual to triangulations. The only change is that 
now the lattice link factor is the two-dimensional scalar propagator, 

f d 2 p eipixi-xj) l 
G(x i ,x j ) = j -^—p-— = —K (rfi\xi-£ j \), (4) 

where Kq is a modified Bessel function. Thus, at the leading order in N, we obtain a sum 
over the planar triangulated random surfaces embedded in two dimensions* The logarithmic 
divergence of G at small separations is very mild. If the tadpole graphs are discarded, as 
they should be because they do not correspond to good triangulations, then each separate 
Feynman graph is finite. This is similar to what we find in the matrix models for c < 1. 
For these types of theories there are general arguments [15] indicating that, for sufficiently 
small A//i, the sum over planar graphs is finite? Therefore, we may look for singularities in 
various physical quantities at some critical value of the dimensionless parameter A///. The 
crucial question is whether, as for c = 1, a sensible string theory can be defined near this 
critical coupling. Although the answer is almost certainly negative, it is interesting to study 
the nature of the critical point. In these notes we will describe our recent numerical results 
which suggest that the spectrum of the theory is tachyonic, and that the string tension 

* Note that, as for c = 1, the coordinates are specified at the centers of the triangles, 
f This is essentially because the number of planar graphs grows only exponentially with the number of 
vertices, and one can put exponential bounds on the value of each graph. 
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actually diverges at the critical point. We will also speculate on how to find a cure for this 
undesirable behavior. 



3. Light-Cone Quantization 



Although the c = 2 matrix model of eq. (3) is certainly more complex than the matrix 
quantum mechanics, we should still be able to take advantage of the simplifications that 
distinguish two-dimensional field theories. The method to pursue this that seems particularly 
convenient, is to continue x° — * ix° and to carry out light-cone quantization [16] of the 
resulting (1 + l)-dimensional field theory, treating x + = (x° + x 1 )/^ as the time and 
x~ = (x° — x x )j\pi as the spatial variable. From the Minkowski signature action, 

S = J dx+dx- Tr (d+Md-M - i/iM 2 + ^=AM 3 ^ , (5) 

we derive the light-cone components of the total momentum 



P + (x + ) = J dx-Tr(d-M) 2 , 
P~(x + )= [dx-Tr(-fiM 2 ^=M 3 ) 



(6) 



Upon quantization, the commutation relations are imposed at equal x + : 

[M i:j {x-),d-M kl {x-)\ = Ud{x- - x~)S a S jk . (7) 

These commutation relations differ from the canonical by the extra factor of ^, which follows 
from the appropriate constrained quantization. This extra factor ensures, for instance, that 
[P + , M] = d-M. Expanding in Fourier modes*, 



Mij 



DO 

^vb/S( a '^ +)e ^ + "'- ( * +)e """) (8) 

o 



| The symbol f is always understood to have purely quantum meaning and does not act on indices. 
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we find the conventional oscillator algebra 



[ aij (k+), a} k (k+)} = 5(k+ - k+)6 a S jk . (9) 
After substituting eq. (8) into eq. (6) and normal ordering, we obtain 



oo 

P + = J dk + k + a\ ) (k + )a ij (k + ) , 


oo 





oo 



X 



+ 







M£k [a^kf + k+)a lk (k+)a kj (k+) + 4(fc+)4^+K#+ + k+)} 
k l k 2 (k 1 + k 2 ) 

(10) 

(repeated indices are summed over). The normal ordering is equivalent to removing the 
tadpole graphs in the Lagrangian approach. P + commutes with P~ so that they can be 
diagonalized simultaneously. As usual in light-cone quantization, the vacuum |0), which is 
defined by 

aij (k + )\0) = , (11) 

is an eigenstate of the fully interacting light-cone Hamiltonian P~ with eigenvalue zero. This 
is related to the fact that the light-cone momentum k + is positive for all quanta. We will 
find that, for sufficiently large A//x, |0) is not the ground state. 

Our goal is to study the ground state and the low-lying excitations. Due to the global 
SU(N) symmetry, M -> StfMQ, the eigenstates can be classified according to their transfor- 
mation properties. We expect the low- lying states to transform as singlets under the SU (N). 
A general singlet state, carrying light-cone momentum P + , can be written as 



oo p * / b \ 

\^(P+)}=J2 / dkidk 2 ...dk h 5 I J2 ki ~ P+ ) 

6=1 n Vi=l / 



P+ 

Ikulho dh,.S iVfc-P+l 

(12) 

fb(h, k 2 , . . . , k b )N' b / 2 TT[aHh) . . . a^k b )]\0) , 
where we have dropped the superscripts + on k{ for brevity. Due to the cyclic property of 
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the trace, the functions f b trivially obey the cyclic symmetry 



f b (ki, fc 2 , ■ ■ ■ , h) = fb(h, fci, • • • , h-\) 



f b (k 2 ,k 3 ,...,k b , fci) 



(13) 



From the normalization condition 



(*(P'+)|tf(P+)) = 5(P+ - P'+) 



we obtain 



J2b / dhdk 2 ...dk b 5 [Y, k i- p+ ] \fb(h,k 2 ,...,k b )\ 2 = l 




(14) 



The singlet states play a special role in this theory: they have no SU(N) degeneracy 
factors and can be thought of as closed strings. Each aj-(fc) creates a string bit carrying 
longitudinal momentum k, and a state of the form 



has the interpretation of a closed string which is a bound state of b string bits. The function 
fb(ki, k 2 , ■ ■ ■ , k b ) superposes states with different distributions of longitudinal momentum 
along the string. If a large- N model is to correspond to a reasonable string theory, it should 
be energetically favorable for all the oscillator indices to be contracted, so that the quanta 
of the matrix field are bound into strings. In fact, if there are finite energy states with free 
indices, then at sufficiently high energy the composite strings are unstable with respect to 
disintegration into bits. As we will discuss in the conclusion, this disease is likely to occur 
in any model with no local gauge symmetry. 

For now, we restrict ourselves to the study of SU (N) singlets, which are the good string 
states. When acting on such states, the light-cone Hamiltonian P" has an important prop- 
erty: in the limit N — >■ oo, P~ takes single closed string states into single closed string states. 
One can easily check that the terms that convert one closed string into two closed strings 



P+ 




dhdk 2 . . . dk b 5 [J2 ki ~ P + ) fb(k h k 2 , ■ ■ ■ , k b )N- b ' 2 Tr[at (fci) . . . a+(A: & )]|0 ) . 



o 
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(two oscillator traces acting on the vacuum) are of order 1/N. Thus, as expected, the string 
coupling constant is ~ 1/-/V, and sending it to zero allows us to study the spectrum of free 
closed string states. The resulting linear light-cone Schroedinger equation can be expressed 
as a set of integral equations for the functions ft, 

b l 

P~fb(h, k 2 ,..., h) = -Mfa, k 2 ,...,k b )^2 — 

1 i=i Ki 

fci ki 

A / f fb+i(Ki k i - k v fc 2, • ■ • , h) | I d] j^ fb+iih, k' 2 , k 2 - k' 2 , fc 3 , . . . , k b ) | 



V^U ^hk'^h - k[) J ^k 2 k' 2 (k 2 - k' 2 ) 

fb-i(ki + k 2 ,k3,...,k b ) fb-i(h,k2 + h,...,k b ) 

\/k x k 2 {kx + k 2 ) y/k 2 k 3 (k 2 + k 3 ) 

(15) 

The terms proportional to A arise either from replacing two neigboring bits by one or from 
dividing a bit into two neigboring bits. The Lorentz invariance of eq. (15) can be made 
explicit if we introduce the longitudinal momentum fractions Xi = ki/P + , and rewrite the 
equation in terms of the functions 

fb(x h x 2 , ...,x b ) = (P+f- 1 V 2 f b (x 1 P+, X 2 P+, . . . , X b P+) . 

From eq. (15) it is not hard to derive 

b l 

2P + P f b (x 1 ,x 2 , ...,x b )= nf b (x!,x 2 , . . . , x b ) V — 



i=l 

Xl _ X2 



A / [ ^ fb+i(x[,xi - x[,x 2 , . . . ,x b ) | f dx r jb+i {x\,x' 2 ,x 2 - x 2 ,x 3 , ...,x b ) | 



2Vn\J ^X\X' X {X\ - x[) J ' y / X 2 X 2 (x 2 - X 2 ) 

fb-l(xi +X 2 ,X 3 ,...,X b ) f b -l(Xh X2 + X3, . . . , X b ) 
y/xiX 2 (xi +X 2 ) ^ X 2 X 3 (X 2 + X 3 ) 

(16) 

and from eq. (14) we find that the functions f b are normalized according to 

oo * / b \ 

J^fe / dxidx 2 . . .dx b 5 I ^Xi - 1 J \f b {xi,x 2 ,. . . ,x b )\ 2 = 1 . (17) 

6=1 q \i=l / 



2P + P is manifestly Lorentz invariant, and its eigenvalues are the squared masses of the 
closed string states. 
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Eqs. (16) and (17) define a rather unusual eigenvalue problem which involves an infinite 
number of unknown functions f b . Each of these functions depends on a number of real 
variables X{ confined to the interval [0, 1], with the constraint Y^=i x i = 1- An important 
issue is the behavior of the functions as X{ — > 0. For /i > the behavior 

lim f b (x 1 ,x 2 ,...,x b ) ^ 

a;i^0 

is inconsistent with a finite energy eigenstate because the contribution ~ ji blows up near 
x\ = faster than the other terms. Guided by this consideration, we will assume that the 
functions f b approach zero as the variables approach the end-points of the interval. 

The system of equations (16) is rather difficult, and we do not know its exact solution. For 
that reason, we will attempt to estimate the spectrum numerically, after introducing a cut- 
off. In light-cone quantization, a particularly convenient method is to replace the continuous 
momentum fractions a; by a discrete set n/K, where the positive integer K is sent to infinity 
as the cut-off is removed [17,18]. Thus, we replace the functions f b (xi,x 2 , ■ ■ ■ ,x b ) by 

g b ( ni , n 2 , ...,n b ) = (K)^/ 2 f b ( ni /K, n 2 /K, n b /K) 



and 



/ 



1 K 



n=i 



According to our discussion of the boundary conditions, we assume that g b is non-vanishing 
only if all its arguments are positive integers (n = is excluded). The cut-off eigenvalue 
problem assumes the form 

9b(ni,n 2 , ...,n b )= Kg b (m,n 2 , . . . , n b ) V — 

is /^SwK.ni -n' 1: n 2 ,...,n b ) ^ g b+ i(n 1: n 2 , n 2 - n' 2 , n 3 , . . . , n b ) 

- M£ + £ — v™>*-«» — + 

9b-i(ni +n 2 ,n 3 ,...,n b ) fl,_i(rai, n 2 + n 3 , . . . , n b ) 
^Jn\n 2 (n\ + n 2 ) \fn 2 n^{n 2 + 77,3) 



(18) 
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where y = A/(2y / 7r/x) is the dimensionless coupling. The normalization condition becomes 



oo K-1 K-1 6-1 



X^SZ'-- \9b(ni,n2,---,n b -i,K -^2rii)\ 2 = 1 . 



(19) 




The cut-off eigenvalue problem reduces to matrix diagonalization because the g's contain 
only a finite number of degrees of freedom, which is equal to the number of partitions of K 
into positive integers, modulo cyclic permutations. All such partitions can be constructed 
recursively with the help of a computer. We find that for K = 10, 11, 12, 13, 14, 15, 16, 17, 18 
the number of degrees of freedom is 107, 187, 351, 631, 1181, 2191, 4115, 7711, 14601 respec- 
tively. As could have been expected, it grows roughly exponentially with K. The maximum 
number of bits is equal to K. As the cut-off K is taken to infinity, the possible values of the 
longitudinal momentum fractions X{ densely populate the interval (0, 1), which corresponds 
to the continuum limit, i.e. the eigenvalue problem of eq. (16). 

If we regard all the independent components of gb(ni, n2, ■ ■ ■ , ^&), with b = 1, 2, . . . , K, as 
components of a vector, then the eigenvalue problem of eq. (18) is equivalent to diagonalizing 
the matrix 



The matrix V is diagonal, while T is off-diagonal, connecting string states whole length 
(number of bits) differs by 1. The entries of V and T can be easily read off from eqs. (18) 
and (19). Alternatively, as shown in ref. 10, V, T, and a basis of normalized states can be 
constructed in terms of creation and annihilation operators. 



The objective of our numerical approach is to compute the low-lying spectrum of the 
matrix (20). First we fix the parameter y, and vary the cut-off K from 9 to 15, which is 
the highest value accessible to us at present. Then, for each y, an extrapolation of the data 
to K — > oo gives us an estimate of the spectrum of the continuum eigenvalue problem of 
eq. (16). Finally we plot the estimated continuum eigenvalues versus y looking for critical 
behavior at some y c . 



2P+P- 



K(V-yT) . 



(20) 



4. Numerical Work 
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What kind of singularity can we expect? Based on our experience with the c = 1 
matrix model, we might expect that the spectrum becomes dense at y — y c , indicating the 
appearance of continuous Liouville momentum p^. The resulting string theory would be 
3-dimensional, with its spectrum given by 

2P + P~ 1 

__^= p 2 +4r __ (21) 

where r runs over non-negative integers, and the tachyonic ground state energy is found 
from the standard formula (2 — D)/6. T is the string tension in the embedding dimensions, 
and pi can be thought of as the "bare" string tension. It is possible that T exhibits some 
dependence on y, and it can even become singular as y — > y c . Depending on the behavior of 
T near the critical point, there are two very different possibilities: 

1) The string tension T is finite at the critical point y c . Then eq. (21) suggests that 
the spectrum of eq. (16) should become continuous at y — y c , starting at a tachyonic value. 
This is the possibility suggested in ref. 10. 

2) T diverges at the critical point, so that all the low-lying eigenvalues of eq. (16) tend 
to — oo. This possibility was not discussed in ref. 10, but the divergence of T has in fact 
been advocated in earlier literature [8]. 

Now we present our current numerical results. Unfortunately, the nature of the numerical 
problem makes it hard to decide with certainty whether 1) or 2) is correct, unless values of 
K considerably higher than 15 can be reached. This is impossible with our present means, 
but may be feasible on a supercomputer. Judging by the available data, 2) is the more likely 
possibility. 

In fig. 1 we plot the lowest three eigenvalues of eq. (18) versus K for y = 0.2. Here 
the convergence to the continuum limit is found to be quite fast. In general, we expect the 
dependence on K to be of the form 



m l (y, K ) = m i(y> °°) + °iM K ~ 



n=l 



In order to estimate the continuum limit, we fitted the dependence on K at fixed y to a ratio 
of two polynomials in 1/K. This fit allows us to extrapolate fig. 1 to arbitrarily large values 
of K. This extrapolation saturates rapidly, showing intuitively expected behavior. 
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In fig. 2 we plot the lowest eigenvalue, m\j versus K for y = 0.2, 0.52, 0.65. These 
plots show the dependence of the rate of convergence on y. For y = 0.52 the convergence is 
much slower than for y = 0.2. For y = 0.65, instead of being concave, the graph becomes 
convex, so that the points no longer appear to converge to a finite limit as K — > oo. 

In fig. 3 we show extrapolations of the lowest eigenvalue to large K. For y = 0.50 (fig. 
3a) the extrapolation converges to a value which is far below that found for K = 15. In fact, 
the continuum value is tachyonic, as anticipated from a D = 3 string theory. For y = 0.55 
(fig. 3b) the extrapolation fails to converge. It follows that there should be a critical value of 
y, which lies between 0.50 and 0.55, where the extrapolated value of m\/ n begins to diverge 
towards — oo. This is indeed what happens. 

In fig. 4 we plot the lowest three eigenvalues, extrapolated to infinite K, versus y. While 
for y = the second and third eigenvalues are degenerate, for y > they are separated by 
a finite gap, which is at first too small to be visible on the figure. The discreteness of the 
spectrum for y > indicates that the quanta are indeed bound into strings. The lowest 
eigenvalue, m\/ )i, begins to dip rapidly for y beyond 0.45. First it becomes tachyonic, and 
then diverges as y approaches 0.53. If we identify this value of y with y c , then the divergence 
of the lowest eigenvalue seems consistent with the divergence of the string tension suggested 
in the possibility 2). Of course, it may be that the gaps in the spectrum, m\ — m\, m\ — m|, 
etc. vanish before the eigenvalues themselves diverge, but fig. 4 indicates otherwise. We find 
that the extrapolated m\f A* plotted vs. y exhibits behavior similar to m\j fi, but diverges at a 
higher value of y, and the gap m\ — ml does not vanish anywhere. If possibility 2) is correct, 
then all the eigenvalues should diverge at the same critical value of y. The fact that this is 
not the case in fig. 4 may be attributed to the lack of precision of our extrapolation. We 
find, however, that the places where m\ and m\ diverge become closer as the extrapolation 
includes the data for higher K. In general, we have to add a cautionary note that the 
precise shape of fig. 4 is not reliable. As y approaches y c , the rate of convergence towards 
the continuum limit becomes slower and slower, so that higher values of K are needed 
for a reliable estimate. Further numerical work is crucial for testing our conjectures. We 
hope, however, that our numerical studies provide a clue about the universal features of the 
continuum limit. 



5. Discussion 
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The possibility 2) suggested in the previous section and supported, to some extent, by 
the available data can be interpreted as an infinite multiplicative renormalization of the 
string tension. In other words, the effective string tension at the critical point y = y c is 
infinitely larger than the bare string tension \i. A very similar scenario was advocated in 
ref. 8 and interpreted as the branched polymer phase of random surfaces. Our results are 
different, though, in that the spectrum at the critical point is tachyonic, as expected from 
the continuum reasoning in eq. (21). 

The divergence of the string tension at y = y c could probably be cured by appropriately 
scaling the bare string tension fi to zero. However, the resulting string theory with finite 
spectrum would still be tachyonic. 

Another problem, which is perhaps even more severe, comes from the states that trans- 
form non-trivially under SU(N), such as the adjoint representation states of the form 

oo p „ + / b \ 

dh... dk b 5 [J2 hi ~ P + ) ft(h, h)N^ /2 l a\ j (k 1 )a] k (h 2 ) ■ ■ ■ aUh^Jk^O) . 

6=1 o \«=1 / 

(22) 

As discussed previously, the SU(N) non-singlet states cannot be identified with closed 
strings. One indication of that is their diverging degeneracy factors in the N — > oo limit. 
Perhaps the non-singlets can be thought of as closed strings that have disintegrated into 
separate bits. Of course, one way of preventing this is through a confinement phenomenon 
which would push their energies to infinity in the continuum limit. In fact, confinement is 
at work in the c = 1 string theory, where the non-singlet energies diverge logarithmically in 
the cut-off [13, 14]. 

Could the c = 2 model of eq. (3) also expel the non-singlet states to infinite energy? The 
answer is negative, as can be seen from a simple variational calculation. Taking a\j{P + )\0) 
as the variational state, and evaluating the matrix element of 2P + P~ / //, gives an upper 
bound of 1 for the lowest eigenvalue. Therefore, there is no confinement, and the non- 
singlet states probably cause additional problems for the interacting string theory. Some 
heuristic arguments with similar conclusions were made in the context of c > 1 models with 
discrete target spaces [13]. The proliferation of the non-singlets was associated with the 
Kosterlitz-Thouless vortices which wind around the target space plaquettes of lattice size. 
An advantage of the light-cone approach is that the non-singlet spectrum can be found with 
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the same technique as the singlet spectrum. It would be quite useful to carry out a more 
detailed numerical study of the spectrum in the adjoint representation. 

The lack of confinement of the non-singlet states can be traced to the fact that the action 
(3) has only the global SU (N) symmetry, which is not gauged. Gauging the SU (N) leads 
to the confinement, and the non-singlets are pushed to infinite energy [19]. The effect of 
gauging the SU (N) symmetry on the two-dimensional theory of eq. (3) is to give the action 

S ga u g ed = J d 2 xTr (JpF^ + \{d a M + i[A a , M]) 2 + i/iM 2 - ^M 3 ) . (23) 

Light-cone quantization of this action was investigated in ref. 19. There a linear light- 
cone Schroedinger equation for the N — > oo limit, similar to eq. (15), was derived. At the 
moment we see no major obstacles to studying this equation using the discretized longitudinal 
momentum cut-off. Perhaps, the theory (23) has new types of critical behavior which lead 
to more interesting string theories than those originating from the model (3). 
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1. The lowest three eigenvalues of eq. (18) versus K for y = 0.2. 

2. The plots of the lowest eigenvalue of eq. (18) versus K for y = 0.2, 0.52, 0.65. 

3. Extrapolation to large K of the lowest eigenvalue for a) y = 0.50; b) y — 0.55. 

4. The lowest three eigenvalues, extrapolated to infinite K, versus y. 
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